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Thermoelectric generators as self-oscillating heat engines 
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In the previous paper [1] a model of a solar cell has been proposed in which the non-periodic 
source of energy - photon flux - drives the collective periodic motion of electrons in a form of plasma 
oscillations. Subsequently, plasma oscillations are rectified by the p-n junction diode into the direct 
current (work). This approach makes a solar cell similar to standard macroscopic heat motors or 
turbines which always contain two heat baths, the working medium and the periodically moving 
piston or rotor. Here, a very similar model is proposed in order to describe the operation principles of 
thermoelectric generators based either on bimetallic or semiconductor p-n junctions. Again plasma 
oscillation corresponds to a piston and sunlight is replaced by a hot bath. The mathematical 
formalism is based on the Markovian master equations which can be derived in a rigorous way from 
the underlying Hamiltonian models and are consistent with the laws of thermodynamics. 


INTRODUCTION 

While the macroscopic heat motors or turbines con¬ 
tain always periodically moving elements like pistons, 
flywheels or rotors [2],[3], the photovoltaic devices, ther¬ 
moelectric generators or biological engines driven either 
by sunlight (photosynthesis) or chemical energy seem, 
according to the standard wisdom, to avoid any self¬ 
oscillation mechanisms. On the other hand extensive 
theoretical studies of quantum heat engines within the 
formalism of quantum open systems [4] -[8] show that the 
process of work extraction must involve a work reservoir 
being a system of a single degree of freedom and execut¬ 
ing an oscillatory motion which in the semi-classical limit 
can be replaced by the periodic external driving. The 
presence of such a system is even necessary to properly 
define the notion of work as a deterministic form of en¬ 
ergy. It suggests that for the examples mentioned above a 
certain “hidden” self-oscillation mechanism of work gen¬ 
eration must also be present. Indeed, as already noticed 
the standard picture of current generation as caused by 
the emerging electric field in a junction cannot be cor¬ 
rect because a DC current cannot be driven in a closed 
circuit by a purely electrical potential difference [9]. In 
the previous paper [1] it was shown that the consistent 
description of work generation by a solar cell can be for¬ 
mulated introducing a work reservoir in a form of THz 
plasma oscillations. 

The aim of the present contribution is to show that 
a very similar model can explain the operation princi¬ 
ples of thermoelectric generators. The basic difference 
between semiconductor and bimetallic devices concerns 
the frequency of plasma oscillations. For doped semicon¬ 
ductors plasma frequencies are in the THz domain (low 
infrared) [10], [11] and can be treated as slow while for 
metals they correspond to the ultraviolet region (PHz 
= 10 15 Hz) and therefore are too fast to be effectively 
coupled to the thermal motion of a lattice. However, it 
is argued in the Appendix, using hydrodynamical model 
of electron gas, that under reasonable assumptions there 
should exist THz modes of plasma oscillations in bimetal¬ 


lic contacts. 


HEAT ENGINE WITH A SLOW PISTON 

Consider a simplified version of the generic quantum 
heat engine model introduced in [4] and successfully used 
in [1] in the context of solar cells. The system, corre¬ 
sponding to the “working medium ” interacts weakly with 
two heat baths at different temperatures. The macro¬ 
scopic piston is described by the external driving V(t) 
added to the free Hamiltonian Hq of the system. 

For slowly varying V (t) (in comparison to the fast in¬ 
ternal dynamics ) the irreversible evolution of the sys¬ 
tem’s time-dependent reduced density matrix p(t ) satis¬ 
fies the following Markovian master equation (MME) 

^Lp{t) = —i[H(t), p(t)\ + C{t)p{t), (1) 

where H(t) = Hq + V(t) is the total Hamiltonian of the 
system and C(t) describes the action of the baths on the 
system ( h = ks = 1 ) • 

The periodic driving executed by a large semi-classical 
oscillator commutes with the system Hamiltonian and 
reads 

V(t) = g(amSU)M, [H o ,M] = 0, (2) 

where j <C 1 is a small amplitude of oscillations. The 
dissipative generator C(t) obtained by the weak coupling 
limit procedure is a function of the magnitude of pertur¬ 
bation and can be written as 

£(*) = £[£(t)], £(*) = g sin fit. (3) 

For all £ the generators £[£] commute and possess the 
stationary states p[£] satisfying the identities 

£[«] = o, = (4) 

where £'[£] = j|£[£] , p'[£] = ^p[£]. 
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Formula for power 


The power P(t) provided by the engine and the net 
heat current J (t) supplied by the baths are defined as [4] 


F(() = -■&(„(*)«), J(t) = TV(5) 

Those definitions are the only ones which are consistent 
with the first and second law of thermodynamics and 
the time-dependence of the Hamiltonian is necessary to 
define work. The stationary average power output of the 
engine reads 


rto 


P = — gQ lim — / 

i 0 ->oo f 0 J 0 


Tr 


(p(f)M) 


cos fit dt. (6) 


Expanding the average power with respect to g one ob¬ 
tains the second order approximation [1] 


P = -g 2 Tr 

2 y 


(p'[o] 


n 2 


S2 2 + £*“[0] 


‘[0 ] M ), 


( 7 ) 


where C* [0] is the Heisenberg picture counterpart of the 
Schroedinger picture generator £[0]. To derive the final 
expression one assumes that, although the modulation is 
slow with respect to the intrinsic motion of the system, 
its frequency is much higher than the relaxation rate of 
the observable M hidden in £*[0]. Then, we have 

P = -\g 2 Tr(p'[0}C*l0]M). (8) 

The lowest order formula (8) is still consistent with ther¬ 
modynamics [1] and is the basic one for the further anal¬ 
ysis of thermoelectric generators. Notice, that the sta¬ 
tionary output power is proportional to the square of the 
amplitude of piston oscillations. This amplitude is a free 
parameter which is determined by the energy flux from 
the hot bath and the load attached to the oscillator (e.g. 
the resistance of the external electric circuit). 


MODEL OF THERMOELECTRIC GENERATOR 

The model of thermoelectric generator is illustrated on 
Fig.l. The working medium is an electron gas distributed 
in two boxes A and B , corresponding to different mate¬ 
rials, metals or doped semiconductors, connected by a 
junction. The electrostatic potential jump through the 
junction denoted by E g is typically of the order of few 
eV. We assume that in the case of metals the density 
of electrons in the box A is lower than in B what im¬ 
plies that potential energy of the electron in the box A 
is higher than in the box B. For semiconductors the box 
A is a p-type and the box B a n-type doped material. 
The hot bath interacts with the electrons in the interface 
region influencing their transport through the junction 



FIG. 1. Model of thermoelectric generator. Bimetallic 
or p-n junction is heated by a hot bath increasing the tran¬ 
sition rate between both regions. The electrons in bulk are 
cooled down by a cold bath. Rectified plasma oscillations pro¬ 
duce a net flow of electrons from the region with their higher 
density (box B) to the lower one (box A). 


while the cold bath at the ambient temperature cools 
down the electrons in the bulk. In the junction region an 
interface between two different material is formed with 
the local concentration of charges and build-in electro¬ 
static potential. One assumes the existence of collective 
charge oscillations localized at the junction with typical 
frequencies in THz domain. For semiconductor devices 
such oscillations are experimentally confirmed and inter¬ 
preted as plasma oscillations with the plasma frequency 


CJp — 



(9) 


where n is a density of charge carriers and m» their effec¬ 
tive mass. The bimetallic junction is more complicated 
as the bulk and surface plasma oscillations in metals pos¬ 
sess frequencies in the PHz domain (10 15 Hz) and hence 
cannot be effectively excited by thermal phonons. In the 
Appendix the mechanism of slow coherent charge oscilla¬ 
tions is proposed which is related to the particular phys¬ 
ical structure of metal-metal junctions. 

The electrons in two boxes A and B are described by 
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two sets of the annihilation and creation operators Ok, 
a k and b\, bj, respectively, subject to canonical anticom¬ 
mutation relations. The electrons are treated as non¬ 
interacting fermions moving in a self-consistent potential 
with the unperturbed Hamiltonian 

Hq = ^ ( E a (kja^dk T E E b (l)b\b h (10) 

k 1 

In a junction a non-homogeneous free carrier distribu¬ 
tion created in a self-consistent build-in potential can be 
perturbed producing collective plasma oscillations with 
the frequency fi. Those oscillations modulate periodi¬ 
cally the Hamiltonian (10) but the detailed mechanism 
is slightly different for the semiconductor and bimetal¬ 
lic junctions (see Appendix for the derivations). In both 
cases the associated time-dependent perturbation which 
should be added to the electronic Hamiltonian (10), is 
proportional to the difference of electron densities in the 
boxes B and A and has a mean-field form (£(f) = g sin f It) 

= HE g Vj (E £ 6+6, - E £ at a k ) . (11) 

1 k 

Here, £ is a small dimensionless parameter describing the 
magnitude of deformation, Va, Vb volumes of the boxes, 
Vj is the effective volume of the junction region and E g is 
the relevant energy scale chosen here to be the potential 
jump across the junction. Remember, that for the semi¬ 
conductor p-n junction the box A corresponds to the “p” 
part and B to the “n” part. In this case —a k a k should be 
rather interpreted as the number operator of holes minus 
an irrelevant constant. 

The interaction with two independent heat bath is de¬ 
scribed by the sum of two terms: C c for the cold bath at 
the ambient temperature T and Ch for the hot bath at the 
temperature Tj > T. The coupling to the cold bath ther- 
malizes electrons in both boxes independently and hence 
preserves separately both electron numbers )T) k a k a k and 
)T k b' k bu ■ It means that C*M = 0 and £* disappears 
from the formula (8) for the stationary power output. 
Notice, that in the discussed idealized model heat trans¬ 
port leading to a temperature gradient is neglected what 
corresponds to the infinite heat conductivity of the bulk. 


Hot bath generator and stationary state 


where 

C$P = 7ki(a k &i P&ifflk ~ 7j{fri a k a k^ P}) 

+ e - ' B »/ Tl 7ki(a] t &i/9&Iak - ^{bjaka^&i, p}), (13) 

The single term (13) describes the following processes: 


a) electron transfer from the box A to B accompanied 
by a positive energy release to the hot bath equal to 


K(k) - E b ( 1) ~ Eg, (14) 


b) the inverse process of electron transfer form the 
lower energy state in B to the higher energy one in A 
with the probability suppressed by the Boltzmann fac¬ 
tor. 

The condition (14) follows from the approximate en¬ 
ergy conservation valid for the tunneling process through 
a barrier of the height E g assisted by thermal fluctuations 
of the hot bath. It means also that the relaxation rates 
7ki are essentially different from zero only if (14) is sat¬ 
isfied. They are proportional to the spectral density of 
the hot bath at w k i = E a ( k) — E b ( 1) — E g . 

Although the generator C* is absent in the formula for 
power the thermal relaxation of electrons in a given box is 
the fastest process and hence determines the form of the 
stationary state p. Within a reasonable approximation 
one can assume that the stationary state of the electronic 
systems with the total Hamiltonian Hq+^M is a product 
of grand canonical ensembles for electrons in both boxes 
with the same temperature T of the device and different 
electro-chemical potentials p a and p b , respectively. The 
associated density matrix has form 


^ ]= zju exp r?E 


E a ( k) + f;E g E- - p a ) aj^k 


+ f Eb (k) - {Eg^-- p b 



(15) 


The electro-chemical potentials p a , Pb are determined by 
the numbers of carriers and hence by densities of elec¬ 
trons in both boxes. The difference of electro-chemical 
potentials is related to the measured voltage $ between 
A and B 


Pa- Pb = 


(16) 


The generator Ch describes thermally induced electron 
transitions from one box to another across the potential 
barrier and can be written in the following form 

^[0]=E £ w (12) 

kl 


Power and efficiency 

One can insert all elements computed in the previous 
section into the expression for power (8). Then one uses 
the properties of the quasi-free (fermionic Gaussian) sta¬ 
tionary state (15) which allow to reduce the averages of 
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even products of annihilation and creation fermionic op¬ 
erators into sums of products of the only non-vanishing 
two-point correlations 

(^k^k^O — ^kk' fc (k), (ukUk/)o — ^kk' (1 ,/c(k) ) , 

(blb^o = Swf 6(1), <6i&J,)o = 6w (1 - f b {l)) • 

Here (• • • )o denotes the quantum average with respect to 
the state p[0] given by (15), and f a { k) and f b ( 1) are the 
Fermi-Dirac statistical distribution functions 

= e 0(E a ( k)- Ma ) + !’ = e ^(B b (l)-Mb) + i’ 

( 18 ) 

with P = 1/T. Inserting the elements defined above into 
the general formula (8) one obtains the leading order con¬ 
tribution to power in a following form 

_ g 2 E 2 V? 

P= £ (V A + V B )(n b -n a )x (19) 

E 7* ( e ~ EJTl t 1 - AO) - [! - /fc(l)] /»(k)), 

kl 


where n a = -^(Sk a k a k)o, n b = y^{J2i b \b i)o are elec¬ 
tron densities. Using (16) one can rewrite (19) as 

~2^>2y2 

P= £ J (V A + V B )(n b -n a )Tx (20) 



where T = Eki7ki[l - /b(l)]/a(k) > 0 is finite in 

the thermodynamical limit. 

The condition for work generation in the discussed 
model of idealized thermoelectric device reads (AT = 
T\ — T) 


e<f> < e<F 0 = E g 




( 21 ) 


The inequality in (21) implies that <f>o can be interpreted 
as an open-circuit voltage of the device and as the 
relative Seebeck coefficient. 

By closing the external circuit one reduces the voltage 
and hence the output power is strictly positive driving 
collective charge oscillations (positive feedback). The 
charge oscillations are subsequently rectified by the diode 
mechanism of the junction producing a direct current. 

The presence of the Carnot factor 1 — T/Ti suggests 
also the interpretation of the eq. (21) in terms of ther¬ 
modynamical efficiency. Indeed, the transport of a single 
electron from the box B to A through the junction re¬ 
quires at least E g of thermal energy extracted from the 
hot bath. Then, a part of energy E g is transformed into 
useful work, equal at most e$o per single electron flowing 
in the external circuit. 

For real systems the efficiency is much smaller than the 
Carnot bound because of the finite heat conductivity and 


damping of plasma oscillations. Similarly, for a typical 
Ti ~ 500/\ , the relative Seebeck coefficient given by (21) 
is of the order of mVjK what is comparable to the high¬ 
est values obtained for some semiconductors. For metals, 
the neglected irreversible transport processes reduce the 
Seebeck coefficient by two or three orders of magnitude. 


CONCLUSIONS 


The aim of this paper is to show that the idea of 
“hidden self-oscillations” proposed in [1] for solar cells 
can be extended to thermoelectric generators. It is illus¬ 
trated by the idealized model in which heat conduction 
is neglected. The “moving piston” is again represented 
by collective charge oscillations with frequencies in the 
THz domain. While such oscillations are observed in 
p-n junctions, for bimetallic ones a plausible mechanism 
of their generation is proposed in the Appendix. The 
challenging open question is to apply the similar ideas to 
systems based on organic molecules including “biological 
engines”. 
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Appendix 

In the Appendix the arguments for the existence of 
slow plasma oscillation modes in bimetallic junctions are 
presented and the origin of diagonal modulation Hamilto¬ 
nian is explained for both, semiconductor and bimetallic 
devices. Here, h is put explicitly in the formulas. 


Slow charge oscillations in bimetallic junctions 


In order to find collective oscillating modes of elec¬ 
tron gas in a bimetallic junction one can use the one¬ 
dimensional hydrodynamical model. The basic variables 
are: the electron density n(x,t), the electrical potential 
tj}{x,t) and the velocity field v(x,t) satisfying the set of 
coupled equations 


dt 


-u—u 
ox 


Id d 

~P ~ e ~' l P = 0 


dx 


dx 


( 22 ) 


d 9 . . „ 

-n+-(n,)= 0 


(23) 



e 

— [n - n 0 ) = 0. 
eo 


(24) 


Taking the derivative of eq.(26) with respect to x 
then inserting relations (27),(28) (omitting the deriva¬ 
tives of n) one obtains the following one-dimensional 
Klein-Gordon-like equation 


■cP_ _ 2 _d 2 _ 

.W CF dx* 



5n = 0 


(29) 


for the density perturbation Sn. Here, the maximal ve¬ 
locity cf is given by 


2 = 5 (3?r 2 ) 2/3 fi 2 2/3 
F 3 2m* 2 



(30) 


where vf is the Fermi velocity and uj p is given by (9). 

For an inhomogeneous material like a bimetallic junc¬ 
tion cf and u> p are assumed to be slowly varying functions 
of x. Using the analogy to the relativistic Klein-Gordon 
equation one concludes that the center of the localized 
perturbation Sn moves like a fictitious one-dinrensional 
particle with the “relativistic” Hamiltonian 


H(x,p) = sJc 2 f (x)p 2 + ft 2 w 2 (x). (31) 

The approximate Hamiltonian equations in the “non- 
relativistic regime” (small p) read 




(32) 


Here m* is the effective mass of an electron and no = 
no(x) is the time-independent density of background pos¬ 
itive ions (jelium model). The pressure p(x 7 t) is not an 
independent variable but is given by the standard for¬ 
mula for the degenerated electron gas 


p = 


(3?r 2 ) 2 / 3 ^ 2 
2m* 


n 5 / 3 = cm 5 / 3 . 


(25) 


The stationary solution is given by v = 0 and the den¬ 
sity fi(x) which up to a certain smoothing essentially fol¬ 
lows the background density no(x). Consider a small 
perturbation of the stationary electron gas distribution 
in a form of well-localized wave packet, i.e. n(x, t ) = 
n(x) + 5n{x,t). Inserting also i(>(x,t) = 4>{x) + Sip(x,t) 
into eqs. (22)-(24) and assuming that the stationary den¬ 
sity is a slowly varying function of x in comparison with 
the variation of the wave packet Sn one obtains a set of 
linearized equations 


m *dt V+ 3nV3 dx 


r\ n 

Sn — e—Si/j = 0, (26) 

dx 


d . _ d 

+ „_„ = o , 


(27) 
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5 ^ - -Sn = 0. 

ax z eo 


( 28 ) 


p=—^U(x), U(x) = Huj p (x). (33) 

To advocate the existence of slow modes in bimetallic 
junctions one can propose the qualitative but plausible 
shape of the effective potential U(x) which interpolates 
between two bulk values characteristic for both metals 
(see Fig. 2). The well in the middle can be explained 
by the presence of a large number of defects in the tran¬ 
sition region between two different lattices. Those de¬ 
fects trap a certain number of electrons reducing the 
density of free electrons and hence also the local value 
of u> p (x) ~ y/n(x). As the jump of the effective potential 
A U = /i(wp B) — w p^) is °f the order of few eV one can 
estimate the depth of the well as A E = 0.1 — leV (still 
higher than the typical value of Tj ~ 0.05eU) and a rea¬ 
sonable value for its width as L = 100 nm (comparable 
to a typical roughness of well polished metallic surfaces). 
Using the harmonic well approximation U(x) ~ |MH 2 i 2 
one can estimate the oscillation frequency of the plas- 
monic wave packet confined in the well by 


n = 


2A E 
ML 2 ' 


(34) 


Taking typical values ~ 10 15 s _1 and vf — 10 62 j one 
obtains M ~ 10 ~ 30 kg and hence, finally Q ~ IT Hz - 
the value comparable to that for plasma oscillations in 
semiconductor p-n junctions. 
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FIG. 2. Bimetallic junction. The shape of the effective 
potential for the plasmonic wave packet. It is instructive to 
compare the proposed shape of U(x), proportional to yjn(x), 
with the very similar shape of electron number density ob¬ 
tained numerically in [12] for the case of Al-Mg contact. In 
[12] (FIG.l) the ideal plane surfaces are separated by a small 
distance of 0.3nm. This gap produces physical effects similar 
to the real metal-metal interface. 


Effective modulation Hamiltonians 


Within this approximation one can use the effective time- 
dependent Hamiltonian with diagonal perturbation 

H(i) = H 0 + (37) 

j 

For the case of bimetallic junction one should take into 
account the fact that the source of perturbation is a 
“wave packet” with fast “internal” charge oscillations at 
the frequency oj p and slow oscillatory motion of its enve¬ 
lope. Hence, the resulting time-dependent Hamiltonian 
possesses the following structure 

H[t) = H 0 + A [/(f) cos u p t] V, (38) 

where f(t) describes the slowly varying envelope. In this 
case the terms linear in A are practically averaged out. 
Using the second order of the standard perturbation the¬ 
ory 


h-tj + x 2^—(0) - TW^k ■ 


(39) 


S # ~£?f ) +A(^ 0) |V r ^ 0) ) + A 2 53 

k^3 


I(4°VI^ 0) )I 2 

Ef ] - 


, (40) 


The aim of this section is to justify the choice of the di¬ 
agonal modulation as a generic consequence of the weak 
coupling between the system and the driving oscillations. 

Consider first the case of slow and weak perturbation 
described by the Hamiltonian Hq + XV(t) which is ap¬ 
plicable to p-n junction. For a fixed t one can apply 
the lowest order of the standard perturbation theory for 
approximate eigenvectors and eigenenergies 

<S>3-<t>f\ (35) 

+ ( 36 ) 


and averaging over fast oscillations one obtains the fol¬ 
lowing effective Hamiltonian 


m = h°+•£ z «f 1 


E (0 ] _ E (0) 


3 k^j 

with, again, diagonal modulation term. 


(41) 


Notice, that for the bimetallic junction the coupling 
of electrons to charge oscillations is the second order ef¬ 
fect and hence, much weaker than for the semiconductor 
devices. However, this is compensated by much higher 
electron densities in metals. 















